Radiation-intensity and temperature dependence of microwave-induced 
magnetoresistance oscillations in high-mobility two-dimensional electron systems 
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We present a detailed theoretical investigation on the radiation induced giant magnetoresistance 
oscillations recently discovered in high-mobility two-dimensional electron gas. Electron interactions 
with impurities, transverse and longitudinal acoustic phonons in GaAs-based heterosystems are 
considered simultaneously. Multiphoton-assisted impurity scatterings are shown to be the primary 
origin of the resistance oscillation. Based on the balance-equation theory developed for magneto- 
transport in Faraday geometry, we are able not only to reproduce the observed period, phase and 
the negative resistivity of the main oscillations, but also to predict the secondary peak/valley struc- 
tures relating to two-photon and three-photon processes. The dependence of the magnetoresistance 
oscillation on microwave intensity, the role of dc bias current and the effect of elevated electron tem- 
perature are discussed. Furthermore, we propose that the temperature-dependence of the resistance 
oscillation stems from the growth of the Landau level broadening due to the enhancement of acoustic 
phonon scattering with increasing lattice temperature. The calculated temperature- variation of the 
oscillation agrees well with experimental observations. 

PACS numbers: 73.50.Jt, 73.40.-c, 78.67.-n, 78.20.Ls 
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I. INTRODUCTION 

Tremendous interest in magneto-transport in two- 
dimensional electron system (2DES) has recently revived 
since the experimental discovery of giant oscillations of 
the longitudinal resistance as a function of the mag- 
netic field in high mohility two-dimensional (2D) elec- 
tron gas (EG) subjected to a microwave radiationfi*^ 
particularly following the recent observations of "zero- 
resistance" states in very clean samples by two indepen- 
dent groupsi 4 i 5 i 6 i 7 These radiation-induced oscillations 
of the longitudinal magnetoresistivity R xx are periodi- 
cal in inverse magnetic field 1/B but, unlike the well- 
known Shubnikov-de Haas (SdH) oscillation whose pe- 
riod is controlled by the electron density N e ^ their pe- 
riods are determined by the radiation frequency u>. The 
observed R xx oscillations exhibit a smooth magnetic-field 
variation with the resistivity maxima at u)/w c = j — <5_ 
and minima at lo/lu c = j + S + (lu c is the cyclotron fre- 
quency, j = 1,2,3...) having positive S± ranging around 
0.1 - 0.25M The resistivity minimum goes downward 
with increasing sample mobility and/or increasing radi- 
ation intensity until a vanishing resistance state shows 
up, while the Hall resistivity keeps the classical form 
R xy = B/eN e with no sign of quantum Hall plateau over 
the whole magnetic field range exhibiting R xx oscillation. 
Later independent experiments^^ confirmed these re- 
sults and the corresponding zero-conductance states were 
also observed in the Corbino samplesiii 

To explore the origin of these peculiar "zero- 
resistance" states, different mechanisms have been 
suggestediiSiiii^iiSiiSiii As is shown by Andreev et aL 13 
that a negative linear conductance implies that the zero 
current state is intrinsically unstable: the system spon- 
taneously develops a non-vanishing local current density 
which is determined by the condition that the compo- 
nent of electric field parallel to the local current vanishes. 



Thus the appearance of negative longitudinal resistivity 
or conductivity in a uniform model suffices to explain the 
observed vanishing resistance. The possibility of absolute 
negative photoconductance in a 2DES subject to a per- 
pendicular magnetic field was first explored 30 years ago 
by RyzhiiiiSiiS Experimentally Keay et aJSireported the 
observation of absolute negative conductance in sequen- 
tial resonant tunneling superlattices driven by intense 
terahertz radiation. Recent worksiSiilii 5 . indicated that 
the periodical structure of the density of states (DOS) 
of the 2DEG in a magnetic field and the photon-excited 
electron scatterings by impurities are the main origin of 
the magnetoresistance oscillations. Durst et. alM> pro- 
posed a microscopic analysis for the conductivity assum- 
ing a 5-correlated disorder and a simple form of the 2D 
electron self-energy oscillatory with the magnetic field, 
obtaining the correct period, phase and the possible neg- 
ative resistivity. Shi and Xie^ gave a similar result us- 
ing the Tien and Gorden current formula™ for photon- 
assisted coherent tunneling. A more quantitative theo- 
retical description was reported recently using a balance 
equation approach developed for radiation-induced mag- 
netotransport in Faraday geometry^ not only reproduc- 
ing the correct period, phase and the negative resistivity 
of the main oscillations, but also predicting the secondary 
peaks and additional maxima and minima observed in the 



experiment 
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and identifying them as arising from 



double- and triple-photon processes. A quantum Boltz- 
mann equation approach based on self-consistent Born 
approximation for large filling factors has also been pre- 
sented very recently, taking account of the elastic (im- 
purity) scattering as the major mechanism for radiation- 
induced magnetoresistance oscillationiS 3 . In addition to 
photon-assisted impurity scattering referred above as the 
mechanism of the absolute negative conductivity (ANC) , 
other possible mechanisms were also explored in the lit- 
erature. Ryzhii et al. proposed that acoustic phonon 
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scatterin gaaaaa and heating of electrons^ could also serve 
as the mechanisms of ANC in 2DES, and made a further 
attempt to connect them with experiments^ 

One of the most important and interesting features of 
this phenomenon is its sensitive temperature dependence. 
The "zero- resistance" states and radiation-induced mag- 
netoresistance oscillations show up strongly only at low 
temperatures typically around T — 1 K or lower. At fixed 
microwave power with increasing temperature, not only 
the zero-resistance regions become narrower and even- 
tually disappear, the whole oscillatory structure (peaks 
and valleys) diminish as well. At temperature T > 4-5 K, 
oscillatory structure disappears completely and the re- 
sistivivty i? 2 . a .-versus-magnetic field becomes essentially 
structureless^^ The temperature-variation of R xx at 
deepest minima exhibits approximate activated-type be- 
havior exp(Xb/T). However, the activation energies ob- 
served by both groups are very high: up to 10 K and 20 K 
at j = 1 minimum respectively^^ These values are about 
an order of magnitude higher than the microwave photon 
energy (u> ^3-5 K) and Landau- level spacing (uj c <2K). 
Furthermore, different To values observed by the two 
groups indicate that the disappearing speed of the os- 
cillatory structure with increasing temperature is sample 
dependent To my knowledge, there has been no theo- 
retical attempt to explain this temperature dependence, 
except that a conjecture of the formation of an energy gap 
around the Fermi surface is suggested under microwave 
irradiation around the resistance minima. 4 

The recently constructed balance-equation model ; 22 ! 2 ? 
which provides a quantitative and tractable approach to 
radiation-induced magnetotransport in Faraday geome- 
try, enables us not only to analyze the magnetoresis- 
tance oscillation, its dependence on the radiation inten- 
sity, but also to deal with its temperature variation in 
a comprehensive way. We suggest that the temperature 
suppression of the magnetoresistance oscillation in these 
high-mobility 2DES comes mainly from the growth of the 
Landau level broadening due to the rapid enhancement of 
acoustic phonon scatterings with increasing temperature 
at this low temperature range. 

In this paper we will carry out a detailed theoretical in- 
vestigation on the different aspects of radiation-induced 
magnetoresitance oscillation. The paper is organized as 
follows. For the convenience and completeness we present 
the general theoretical model and formulation in Section 
II. As a typical example we analyze GaAs-based systems 
for which the relevant material parameters are discussed 
in the Section III. Section IV concentrates on the trans- 
port properties of GaAs-based heterosystems at lattice 
temperature T = IK. We will give a detailed discussion 
on the impurity and acoustic phonon scattering related 
linear and nonlinear magnetoresistance induced by the 
irradiation of 0.1 THz microwave of different intensities, 
and the effect of elevated electron temperature. Section 
V is devoted to the analyses on the lattice-temperature 
dependence of the magnetoresistance oscillation. Finally, 
a brief summary is given in Section VI. 



II. FORMULATION 

A. Balance equations in crossed electric and 
magnetic fields 

The experiments allow us to assume the 2DEG being 
in extended states over the magnetic field range relevant 
to this phenomenon. For a general treatment, we con- 
sider N e electrons in a unit area of a quasi-2D system in 
the x-y plane with a confining potential V(z) in the z- 
direction. These electrons, besides interacting with each 
other, are scattered by random impurities/disorders and 
by phonons in the lattice. To include possible clliptically 
polarized microwave illumination we assume that a uni- 
form dc electric field Eo and a high-frequency (HF) ac 
field of frequency w, 

E t = E s sin(tctf) + E c cos(atf), (1) 

are applied in the x-y plane, together with a magnetic 
field B = (0,0, B) along the z direction. In terms of the 
2D center-of-mass momentum and coordinate of the elec- 
tron system ) 29 i 30 i 31 which are defined as P = J2j Pj\\ an d 

R = N ^ J2j r j\\ with Pill = (Pjx,Pjy) and r,-j| = (a^y,-) 
being the momentum and coordinate of the jth electron 
in the 2D plane, and the relative electron momentum 
and coordinate p^.,, = Pj|| — P/A c and r^.,, = — R. 
the Hamiltonian of the system can be written as the sum 
of a center-of-mass part H cm and a relative electron part 
H el (A(r) is the vector potential of the B field), 

H cm = -^(P - A c eA(R)) 2 - A e e(E + E t ) • R, (2) 
2N c m 

^ = E[^(^n- eA Kii)) 2 + ^+^)] 

+ '£V c (r' q -r' jll ,z u z j ), (3) 

together with electron-impurity and electron-phonon in- 
teractions H c i and H op . Here m and m z are respec- 
tively the electron effective mass parallel and perpendic- 
ular to the plane, and V c stands for the electron-electron 
Coulomb interaction. It should be noted that the uni- 
form electric field (dc and ac) appears only in H cm , and 
that H er is just the Hamiltonian of a quasi-2D system 
subjected to a magnetic field. The coupling between the 
center-of-mass and the relative electrons exists via the 
electron-impurity and electron-phonon interactions. Our 
treatment starts with the Heisenberg operator equations 
for the rates of changes of the center-of-mass velocity 
V = -i[V,H]+dV/dt, with V = -i[R,H], and of the 
relative electron energy H er = — i[H el ,H], and proceeds 
with the determination of their statistical averages. 

As proposed in Ref. l2^. the center-of-mass coordinate 
R and velocity V can be treated classically, i.e. as 
the time-dependent expectation values of the center-of- 
mass coordinate and velocity, R(t) and V(t), such that 



3 



R(t) - R(i') = f t , V(a)da. We are concerned with the 
steady transport state under an irradiation of single fre- 
quency and focus on the photon-induced dc resistivity 
and the energy absorption of the HF field. These quan- 
tities are directly related to the time-averaged and/or 
base-frequency oscillating components of the center-of- 
mass velocity. Although higher harmonics of the current 
may affect the dc and lower harmonic terms of the drift 
velocity through entering the damping force and energy 
exchange rates in the resulting equations, in an ordinary 
semiconductor the power of even the third harmonic cur- 
rent is rather weak as compared to the fundamental. For 
the HF field intensity in the experiments, the effect of 
higher harmonic current is safely negligible. Hence, it 
suffices to assume that the center-of-mass velocity, i.e. 
the electron drift velocity, consists of a dc part Vo and a 
stationary time-dependent part v(t) of the form 



V(t) = vo + vi cos(wi) + V2 sin(wt). 



(4) 



With this, the exponential factor in the operator equa- 
tions can be expanded in terms of Bessel functions 3 n (x): 

e -iq-J t , V(s)ds _ j2 ^ e i(q.v -na>)(W) + 

n— — OO 

OO 

^e™^ Jn(C)Jn-m(Oe i(q-VD "" a,)( ^ t ' ) - 



n— — 00 



Here the argument in the Bessel functions 



[tan 



■ v i) 2 + (qy • v 2) 



(5) 



and tan ip = (q • V2)/(q • Vi). On the other hand, for 2D 
systems having electron sheet density of order of 10 15 
m~ 2 , the intra-band and inter-band Coulomb interac- 
tions are sufficiently strong that it is adequate to describe 
the relative-electron transport state using a single elec- 
tron temperature T c . Except this, the electron-electron 



interaction is treated only in a mean-field level under ran- 
dom phase approximation (RPA) . 29 3( i For the determina- 
tion of unknown parameter Vo, Vi, V2, and T e , it suffices 
to know the damping force up to the base frequency oscil- 
lating term F(t) = F + F s sin(wt) + F c cos(wt), and the 
energy-related quantities up to the time-average term. 
We finally obtain the force and energy balance equations: 



= 
vi = 

-V2 = 



JV o eE + JV o e(v x B) + F , 
eE. s F, e 

(v 2 x B), 



mui N c moj 
eE c , F c 
muj 



e 



N c muj muj 



(vi x B), 



N c eE Q • v + S p - W = 0. 



(6) 
(7) 

(8) 



(9) 



Here 



F = ]T|/y(q||)| 2 J2 q||J'(Ona(q||,wo-nw) 

qn n— — oo 

11 00 

+ ]T|A/(q)| 2 q||J 2 (OA 2 (q^o + ^ q -^),(10) 

q n= — 00 

is the time-averaged damping force, is the time- 
averaged rate of the electron energy-gain from the HF 
field, ^N c e(E s • v 2 + E c • vi), which can be written in a 
form obtained from the right hand side of Eq. IjlUI) by 
replacing the qy factor with nu>, and W is the time- 
averaged rate of the electron energy-loss due to coupling 
with phonons, whose expression can be obtained from the 
second term on the right hand side of Eq. (fTTTf> by replac- 
ing the qy factor with f2 q , the energy of a wavevector- 
q phonon. The oscillating frictional force amplitudes 
F s = F22— F u andF c = F 2 i+Fi 2 are given by (// = 1,2) 



OO OO 

Fi^-EwMqil)! 2 E [Jn(0]'ni(q||,a;o-n W )-^q||^|M(q)| 2 £ [#($] ' Axfa, wo + O q - nw), (11) 

qii n— — oo q n— — oo 

00 00 

F 2M = ]Tq||^|C/(q||)| 2 2nJ 2 (e)H 2 (q||,c -nc) + ^q |t ^|Af( q )| 2 ]T 2nJ 2 1 (£)A 2 (q, u Q + tt q - ulo). (12) 

qj| n— — 00 q n— — oo 

1 



In these expressions, 77^ = qn ■ v^/uj^; ujq ee qn • 
vo; C/(q||) and M(q) stand for effective impurity and 
phonon scattering potentials, n 2 (q||,f2) and A 2 (q, f2) = 
2n 2 (q|| , fi)[n(fi q /T) - n(fl/T e )} (with n(x) = l/(e x - 1)) 
are the imaginary parts of the electron density correlation 
function and electron-phonon correlation function in the 



presence of the magnetic field. Hi (qn , f2) and Ai(q, fi) 
are the real parts of these two correlation functions. 

The HF field enters through the argument £ of the 
Bessel functions in F , F MJ/ , W and S p . Compared with 
that without the HF field (n — term only)^ we see 
that in an electron gas having impurity and/or phonon 
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scattering (otherwise homogeneous), a HF field of fre- 
quency lu opens additional channels for electron transi- 
tion: an electron in a state can absorb or emit one or 
several photons and scattered to a diffcrnt state with 
the help of impurities and/or phonons. The sum over 
\n\ > 1 represents contributions of single and multiple 
photon processes of frequency-o; photons. These photon- 
assisted scatterings help to transfer energy from the HF 
field to the electron system (S p ) and give rise to addi- 
tional damping force on the moving electrons. 

Eqs. form a closed set of equations for the de- 

termination of parameters vq, vi V2 and T c when Eo, E c 
and E s are given in a 2D system subjected to a magnetic 
field B at temperature T. Thus they provide a com- 
prehensive and quantitative description of transport and 
optical properties of magnetically-biased quasi-2D semi- 
conductors subjected to a dc bias and a HF radiation 
field in Faraday geometry. 

Note that vi and V2 always exhibit cyclotroresonance 
in the range ui ~ lu c , as can be seen from Eqs. {71) and © 
rewritten in the form 



ately get the transverse and longitudinal resistivities 



vi 



{- 

I mtu 



(lu 2 — lu 2 ) I mtu 
1 



N e muj 

2 



LU 



v 2 = 



{- 



(lu 2 — lu 2 ) I mtu 
1 



E s + (E c xB) 

mtu 



— (F c xB) 



E c (E s x B) 

TYILU 



N r mLU - 



-(F s x B) 



(13) 



(14) 



Therefore, the argument £ may be significantly different 
from that of the corresponding Bessel functions in the 
case without a magnetic field or with a magnetic field 
in Voigt configuration, where the electron motion is not 
affected by the magnetic fields On the other hand, im- 
purity and phonon scatterings can affect £ through the 
damping forces F s and F c . Eqs. (|13fl and (|14fl when ne- 
glecting the damping forces F s = = F c , yield a Bessel- 
function argument £ equivalent to that used in the early 
literatur o 33 i 34 . The approximation of neglecting damp- 
ing forces is valid only in weak scattering limit and away 
from cyclotroresonance. Depending on Vi and v 2 , the 
damping forces F s and F c in Eqs. (|13|) and (|14fl are im- 
portant not only in the general scattering case over the 
whole magnetic-field range but also in weak scattering 
case in the vicinity of cyclotron resonance in that they 
remove the divergence and yield finite oscillation veloci- 
ties Vi and v 2 at lu = lu c . 



B. Longitudinal and transverse resistivities 

The nonlinear resistivity in the presence of a high- 
frequency field is easily obtained from Eq. ijfJJl. Taking 
vq to be in the x direction, vq = (uoa;,0,0), we immedi- 



_ E 0y = B_ 

xv N c ev 0x iV c e' 



Fn 



Eq x 

N c ev 0x N 2 e 2 v 0x 



(15) 
(16) 



We see that in the present model (with parabolic energy 
spectrum) the transverse resistivity R xy remains the clas- 
sical form, with no change in the presence of HF radiation 
and/or dc current. The longitudinal resistivity R XXl on 
the other hand, can be strongly affected by the irradia- 
tion also by dc bias. The linear longitudinal resistivity is 
the weak dc current limit (vq x — > 0): 



Rxx 



lt%ll)l 2 00 

N?e 2 



E J ««) 



dU 2 



n— — oo 
,2 oo 



dfl 
dA 2 



(17) 



Note that although according to Eas. (|1()J) . (|16J) and 
(|17fl . the linear and nonlinear longitudinal magnetoresis- 
tivity R xx can formally written as the sum of contribu- 
tions from various individual scattering mechanisms, all 
the scattering mechanisms have to be taken into account 
simultaneously in solving the momentum- and energy- 
balance equations 0, (jHJl and © for vi, v 2 and T , 
which enter the Bessel functions and other parts in the 
expression of R xx . 



C. Landau- level broadening 

In the present model the effects of interparticle 
Coulomb interactions are included in the electron com- 
plex density correlation function n(q|| , f2) = LTi (q|| , fi) + 
in2(q||,f2), which, in the random phase approximation, 
can be expressed as 



n(q,|,0) 



n Q (qn,n) 
c(qn»n) 



where 



e(qH , O) = 1 - )n (q|| , O) 
is the complex dynamical dielectric function, 

^2 



V(q 



■H(q 



(18) 



(19) 



(20) 



di- 



is the effective Coulomb potential with n the 
electric constant of the material and ^(^n ) a 
2D wavcfunction-related overlapping integration^ 
n (q||,n) = n i(q||,n) + in 2(q||,O) is the complex 
density correlation function of the independent electron 
system in the presence of the magnetic field. With this 
dynamically screened density correlation function the 
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collective plasma modes of the 2DES are incorporated. 
Disregard these collective modes one can just use a 
static screening e(q|| , 0) instead. 

The II02 , ^) function of a 2D system in a magnetic 
field can be written in terms of Landau representation: 31 

n 02 (q||,O) = ^£c n ,„,(i|gf/2)II 2 (n,n',fi), (21) 



n 2 (n, n', O) = — / ds [/(e) - /(e + SI)] 



x ImG„(e + Sl)lmG n '(e) 



where Ib = -y/l/|e_B is the magnetic length, 
C n>n+l (Y) = n\[(n + ly.r'Y'e^iLKY)} 



(22) 



(23) 



with L l n (Y) the associate Laguerre polynomial, f(e) = 
{exp[(e — /i)/T c ] + l} 1 the Fermi distribution function, 
and ImG rl (e) is the imaginary part of the electron Green's 
function, or the DOS, of the Landau level n. The real 
part function IIoi(q||,r2) and corresponding Aoi(q||,f2) 
function can be derived from their imaginary parts via 
the Kramers-Kronig relation. 

In principle, to obtain the Green's function imG n (e), a 
self-consistent calculation has to be carried out from the 
Dyson equation for the self-energy with all the impurity, 
phonon and other scatterings included. The resultant 
G n is generally a complicated function of the magnetic 
field, temperature, and Landau-level index n, also de- 
pendent on the relative strengths of different kinds of 



scatterings 



35.36 



In the present study we do not attempt 



a self-consistent calculation of G„ (e). Instead, we choose 
a Gaussian-type form— for the purpose of demonstrat- 
ing the observed oscillations (e n is the energy of the n-th 
Landau level): 



electron sheet density N c = 3 x 10 11 cm~ 2 but hav- 
ing linear mobility /io(lK) = 2.4 x 10 7 cm 2 V _1 s _1 and 
^o(lK) = 1.46 x 10 7 cm 2 V _1 s _1 respectively in the 
absence of the magnetic field. In GaAs/AlGaAs sys- 
tems, phonon modes and electron-phonon couplings are 
well established^ We consider both transverse acous- 
tic phonons (interacting with electrons via piezoelectric 
coupling) , and longitudinal acoustic phonons (interacting 
with electrons via piezoelectric and deformation potential 
couplings). We assumed that the elastic scatterings are 
either due to the remote charged impurities which are lo- 
cate a distance s = 60 nm away from the interface of the 
heterojunction in the barrier side, or due to background 
impurities which are evenly distributed throughout the 
GaAs region. The impurity densities are determined 
by the requirement that electron total linear mobility 
equals the giving value at temperature T = 1 K. The 
relevant effective impurity scattering potentials |[/(q||)| 2 
and electron-phonon matrix elements |M(q)| 2 were dis- 
cussed earlieriSi The material and coupling parameters 
for the system are well defined and taken as: electron 
effective mass m = 0.068 m (m e is the free electron 
mass), transverse sound speed v st — 2.48 x 10 3 m/s, lon- 
gitudinal sound speed v s \ — 5.29 x 10 3 m/s, acoustic de- 
formation potential 5 = 8.5 eV, piezoelectric constant 
ex4 = 1.41 x 10 9 V/m, dielectric constant k = 12.9, ma- 
terial mass density d = 5.31 g/cm 3 . The depletion layer 
charge number density is taken as JVd op = 5x 10 10 cm~ 2 . 

In GaAs systems at low temperatures Hq(T) comes 
from impurity, transverse and longitudinal acoustic 
phonon scatterings: 



1 

Mo 



a4 1} 



(pt) 

Mo 



(pi) ' 
Mo 



(26) 



imG n (e) = ~\/^2 ex P 



(e - e n ) 2 



2L 2 



with a broadening width giving by 



/ 2euj c a \ 1 / 2 
\7rm/io(T) 



(24) 



(25) 



where fJ-o(T) is the linear mobility at temperature T 
in the absence of the magnetic field and a > 1 is 
a semiempirical parameter to take account the differ- 
ence of the transport scattering time determining the 
mobility /i (T), to which larger angle scattering con- 
tributes a heavier weight, from the single particle life- 
time, to which scattering with small or large angle equally 
contributes 



III. GA AS-BASED 2DES 

To contact with recent experiments, we focus our at- 
tention to two ultra high mobility two-dimensional elec- 
tron systems of GaAs/AlGaAs heterostructure with same 



The longitudinal magnetoresistivity R xx can also be for- 
mally written as the sum of the impurity, transverse and 
longitudinal acoustic phonon contributions: 



Rxx — R^r\. ~\~ Rt—L ^ ~\~ R^ r . 



(27) 



In the following sections we will carried out numerical 
calculations for Rxl, R^x and RxlP assuming linearly 



polarized MW fields (E c 
included. 



0) with multiphoton processes 



The microwave field intensity required for the appear- 
ance of resistivity oscillation in these high-mobility sam- 
ples is moderate. The slight electron heating induced 
by the irradiation in these systems is unimportant as far 
as the main phenomenon is concerned. As can be seen 
later that the radiation-induced magnetoresistance oscil- 
lation is not sensitive to the moderate rise of electron 
temperature. Therefore, R xx can be obtained directly 
from Eq. |T7|) with T c = T. We will check the effect of 
elevated electron temperature in the subsection IV-c. 



IV. R xx AT TEMPERATURE T = 1 K 
A. Impurity-induced magnetoresistivity 

1. Linear resistivity 

We first assume that the elastic scatterings are due to 
ionized remote impurities. 30 Fig. la shows the impurity- 
induced longitudinal resistivity Rxl versus u>/uj c = 7 C 
subjected to a microwave radiation of frequency lo/2-k — 
0.1 THz at four values of amplitude: E s = 20, 45, 65 and 
80V/cm. At temperature T = IK SdH oscillations of 
period j c = 0.039 show up strongly at high uj c side, and 
then gradually decay away as \/ui c increases. In addi- 
tion to this all four resistivity curves exhibit clear oscil- 
lation having the main oscillation period 7 C = 1 (they 
are crossing at integer points j c = 2,3,4,5). The resis- 
tivity maxima locate around 7 C = j — 5_ and minima 
around j c = j + 5 + with 5± ~ 0.22 - 0.26 for j = 3, 4, 5, 
S± ~ 0.17-0.22 for j = 2, and S± ~ 0.08-0.14 for j = 1 
(Fig. lb) . The amplitude of the oscillation increases with 
increasing HF field intensity for j c > 1.5. Resistivity gets 
into negative value for E s = 80 V/cm around the minima 
at j = 1, 2 and 3, for E s = 65 V/cm at j = 1 and 2, and 
for E s — 20 and 45 V/cm at j = 1. These main peak- 
valley structures are related to single-photon (|n| = 1) 
processes. In the vicinity of j c = 1, where the cyclotron 
resonance greatly enhances the effective amplitude of the 
HF field in photon-assisted scatterings, multiphoton pro- 
cesses show up. The amplitudes of the j — 1 maximum 
and minimum no longer monotonically change with field 
intensity. Furthermore, there appears a shoulder around 
7 C = 1.5 on the curves of E s = 45 and 65 V/cm, and 
it develops into a secondary peak in the E s = 80 V/cm 
case. This peak- valley structure around 7 C = 1.5 is re- 
lated to two-photon = 2) processes. The oscillatory 
peak-valley structure related to three-photon processes 
was also demonstrated in the cases of lower microwave 
frequency (to = 60 and 40GHz)i22 

The dependence of the resistivity at maxima and min- 
ima on the microwave intensity is shown in Fig. lc, where 
we plot the calculated photoresistivity, i.e. the magne- 
toresistivity in the presence of radiation, Rxl , minus the 
dark resistivity Rxl{E s — 0), as a function of radiation 
field amplitude E s at peaks and valleys of j = 2 and 
j = 3. We see that \R XX — Rxx(0)\ grows like E 2 S at lower 
intensity and become linearly dependent on E s at higher 
intensity within the amplitude range shown. This is in 
agreement with experiments4i^ 

Fig. Id shows the positive parts of Rxl /Rxl (E s = 0) at 
minima of j = 2,3 and 4 on logarithmic scale as functions 
of El 

It should be noted that the above feature of impurity- 
induced magnetoresistivity is not very sensitive to the 
form of elastic scattering potential. In Fig. 2a we demon- 
strate Rxl due to background impurity scattering as a 
function of u /lj c under several different microwave ampli- 
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FIG. 1: (a) The longitudinal linear magnetoresistivity R^l 
induced by remote-impurity scattering in a GaAs-based het- 
erosystem subjected to crossed magnetic fields B and in- 
plane linearly polarized HF fields E s sm(ui) of frequency 
lu/2-k = 0.1 THz with several different amplitudes at lat- 
tice temperature T = IK. ui c = eB/m stands for the cy- 
clotron frequency. The other parameters are: electron density 
iV c = 3.0 x 10 11 cm~ 2 , zero-magnetic-field linear dc mobility 
/io(lK) = 2.4 x 10 7 cm 2 V -1 s~\ and the broadening coeffi- 
cient a — 12. The electron temperature is set to be T c — T. 
(b) Parameters 5+ and 5- for locations of resistance maxima 
and minima at several HF field amplitudes, (c) The photore- 
sistivity Rxl — Rxl(E s — 0) at maxima and at minima of 
j = 1 and j = 2 against the amplitude of the HF field, (d) 
Rxl/ Rxl(E s = 0) is shown against E 2 S on logarithmic scale 
for j = 2, 3 and 4. 
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FIG. 2: (a) The linear magnetoresistivity Rxl induced by 
background-impurity scattering in a GaAs-based 2DEG sub- 
jected to linearly polarized HF fields E s sin(wt) of frequency 
lu/2-k = 0.1 THz. The system parameters are the same as in- 
dicated in Fig. 1. (b) The background impurity-induced pho- 
toresistivity Rxl — Rxl (0) at maxima and at minima of j — 1 
and j = 2 shown in Fig. 2a is plotted against the amplitude 
of the HF field. 



7 







1 1 1 1 1 1 1 1 1 1 1 ■ ■ 1 1 1 1 1 

A to/2jt=0.1 THz V ] 
r\ E =80 V/cm --- v =.015v F ; 

L\ t=ik v ° = r vp 

\r\ \ v = 05 v 

y \\ o f 

/ \! A v - 1v f 








\J Remote impurity scattering 



1 2 3 4 5 

Co/co 

c 



FIG. 3: Remote-impurity induced nonlinear longitudinal 
magnetoresistivity R xx as defined in Ea. l|16|l in the GaAs- 
based 2DEG subjected to crossed magnetic fields B and an 
in-plane linearly polarized HF field E s sin(a;£) of frequency 
lu/2-k — 0.1 THz and amplitude E s = 80 V/cm under several 
different dc bias velocities vq = 0, .015, .025, .05 and .Ivf, 
where vf ~ 2.4 x 10 5 m/s is the electron Fermi velocity. 
ui c = eB/m is the cyclotron frequency. The other parame- 
ters are the same as indicated in Fig. 1. 

tildes. Although the scattering potential of background 
impurities appears quite different from that of remote 
impurities ipil the main features of the radiation-induced 
magnetoresistivity oscillation are in common for both 
scattering potentials. The difference shows up only in 
some quantitative details. The dependence of the back- 
ground impurity induced resistivity at maxima and min- 
ima on the microwave intensity shown in Fig. 2b, is also 
similar to the case of remote impurity scattering (Fig. lc). 

2. Nonlinear resistivity 

Fig. 3 shows the nonlinear longitudinal resistivity R xx 
due to remote-impurity scattering calculated directly 
from Eq. iJTSJl for the 2DES subject to a 0.1 THz mi- 
crowave radiation of amplitude E s — 80 V/cm under dif- 
ferent dc bias velocities v = 0, .015, .025, .05 and 0.1 vp, 
where the electron Fermi velocity vf = 2.4 x 10 5 m/s. For 
the given strength of the radiation field the linear magne- 
toresistivity (vanishing Eq or vo) exhibits the strongest 
oscillation. A finite dc bias always suppresses the oscilla- 
tion and may destroy the negative resistivity appearing 
at vanishing dc bias. 

B. Acoustic-phonon-induced magnetoresistivity 

Acoustic phonon scattering has recently been proposed 
as a mechanism of the absolute negative resistivity lead- 
ing to vanishing resistance state^S^ To check such a 
possibility we show in Fig. 4 the linear magnetoresis- 
tivity R XX1 contributed separately by transverse acous- 
tic phonon scattering (a) and by longitudinal acous- 
tic phonon scattering (b), as functions of lu/lu c under 
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FIG. 4: Linear magnetoresistivities Rxx (a) and Rxx (b) 
induced by transverse and longitudinal acoustic phonons in a 
GaAs-based 2DES subjected to microwave fields E s sin(u;i) of 
frequency u)/2ir = 0.1 THz having different amplitudes. The 
lattice temperature is T — 1 K and the electron temperature 
T c — T. The other parameters are: electron density iV c = 
3.0 x 10 n cm- 2 , dc mobility W (1K) = 2.4 x 10 7 cm 2 V" 1 
s _1 , and the broadening coefficient a — 12. 



0.1 THz microwave illumination of different strengths. 
Photon-assisted acoustic-phonon scattering itself indeed 
can give rise to pronounced resistance oscillation with 
changing magnetic field and R xx at oscillation minima 
can go down to negative under sufficiently strong mi- 
crowave irradiation for both transverse and longitudinal 
phonon scatterings. These acoustic-phonon induced R xx 
oscillations, though can still be considered nearly peri- 
odical in inverse magnetic field, exhibit quite different 
behavior from that of impurity-induced R xx oscillation 
shown in Fig. 1 and also different from each other, espe- 
cially in respect to the phase of the oscillation. In the case 
of impurity scattering the resistance maxima always ap- 
pear at 7 C = j — <5_ and the maxima at j c = j + 5 + , while 
in the case of longitudinal acoustic scattering the max- 
ima appear at 7c = J + <5+ and the minima at 7 C = j — r5_ 
except the case of j = 1. The phase of transverse acoustic 
phonon induced R xx oscillation resembles that of impu- 
rity at j = 1,2 but becomes that of longitudinal acoustic 
phonon at j = 3, 4, 5. 

As in the case of impurity scattering a finite dc bias 
always suppresses the phonon-induced resistance oscilla- 
tions and may destroy the negative resistivity appearing 
at vanishing dc bias as seen in Fig. 5, where we plot the 
nonlinear longitudinal resistivity R x p x ^ induced by trans- 
verse acoustic phonons and R x p x ^ induced by longitudi- 
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FIG. 5: Nonlinear magneto-resistivities R^ (a) and R x p ^ 
(b) induced by transverse and longitudinal acoustic phonons 
under different dc bias vo in a GaAs-based 2DES exposed 
to microwave fields E s sin(aji) of frequency u}/2tt = 0.1 THz 
having amplitude E a — 170 V/cm and E s = 95 V/cm respec- 
tively. Here v st = 2.48 x 10 3 m/s and u s] = 5.29 x 10 3 m/s 
are the transverse and longitudinal sound speed. The other 
parameters are the same as indicated in Fig. 4. 

nal acoustic phonons, for the 2DEG subject to a 0.1 THz 
microwave radiation with fixed amplitudes under several 
different dc current biases. These results are in qualita- 
tive agreement with a recent report^ 

However, at temperature T = 1 K the acoustic phonon 
scattering contributes a part of R xx which is more than 
an order of magnitude smaller than that contributed from 
impurity scattering in the system having a mobility of 
2.4 x 10 7 cm 2 V- 1 s _:L at T = IK. Therefore, acoustic 
phonon scattering essentially gives no direct contribu- 
tion to the experimentally observed resistance oscillation 
at this low temperature. Nevertheless, acoustic phonons 
play a key role in suppressing the resistance oscillation at 
elevated lattice temperatures. We will discuss this issue 
in Section V. 




FIG. 6: The impurity-induced linear magnetoresistivity R^l 
(a) and transverse acoustic phonon induced linear magne- 
toresistivity of a GaAs-based 2DES subjected to a mi- 
crowave field of amplitude E s = 80 V/cm and frequency 
Lu/2-K = 0.1 THz at elevated electron temperatures. The lat- 
tice temperature is T = IK. The other parameters are the 
same as indicated in Fig. 4. 



lattice temperature T = 1 K subject to a microwave 
irradiation of frequency oj/2tt — 0.1 THz and ampli- 
tude E s = 80kV/cm. We see that at T c = 5K (even 
much lower) SdH oscillations disappear completely, but 
the radiation-induced oscillations remain essentially the 
same. Only when the electron temperature becomes 
much higher, e.g. 50 K, can appear the appreciable 
change in the resistance oscillation curve. This indicates 
that the radiation-induced resistivity oscillation is quite 
insensitive to electron temperature and we can analyze 
the moderate-strength microwave induced resistivity os- 
cillation by neglecting the electron temperature change 
in the system. On the other hand, using a slightly el- 
evated electron temperature provides a way to separate 
the radiation-induced oscillation from the SdH effect. 

Fig. 6b shows the effect of elevated electron temper- 
ature on transverse acoustic phonon induced resistiv- 
ity Rx P x^ ■ With electron temperature increasing from 
T c = 1 K to 5 K, the magnitude of the oscillation some- 
what grows. But its absolute contribution to magnetore- 
sistivity is still an order of magnitude smaller than the 
of impurity as long as the lattice temperature remains 
T = IK. 



C. Effect of elevated electron temperature 

Upon the microwave irradiation the electron temper- 
ature can be higher than the lattice temperature. To 
have an idea of how the elevated electron temperature af- 
fects the magnetoresistivity oscillation we show in Fig. 6a 
the background impurity induced longitudinal resisitiv- 

ity Rxx of the 2DEG having several different electron 
temperatures T c = 1,5 and 50 K but all at the same 



Experiments showed that the strongest SdH oscilla- 
tions always show up on the dark resistivity curve and 
with enhancing the radiation intensity the SdH oscilla- 
tions weakenAii This fact can be easily understood as 
due to the rise of the electron temperature caused by the 
microwave illumination: SdH oscillations are suppressed 
by the rising electron temperature, while the radiation- 
induced R xx oscillations remain the same as long as the 
lattice temperature keeps unchange. 



V. EFFECT OF ELEVATED LATTICE 
TEMPERATURE 

One of the most important aspect of the exper- 
imental finding on the magnetoresistance oscillation 
in irradiated 2DES is the temperature-dependence of 
R xx . Experiments^ found that the "zero-resistance" 
states and radiation-induced magnetoresistance oscilla- 
tions show up strongly only at low temperatures typ- 
ically around T = 1 K or lower. At fixed microwave 
power with increasing temperature, not only the zero- 
resistance regions become narrower and eventually disap- 
pear, the whole oscillatory structure (peaks and valleys) 
diminish as well. At temperature T >4-5K, oscillatory 
structure disappears completely and the resistivivty R xx - 
versus magnetic field becomes essentially flat£ 

Both groups analyzed the temperature-variation of 
R xx at deepest minima using an activated-type depen- 
dence exp(To/T). The activation energies observed by 
both groups are very high and different: up to 10 K 
and 20 K at j = 1 minimum respectively^^ These val- 
ues are about an order of magnitude higher than the 
microwave photon energy (uj ^3-5 K) and Landau-level 
spacing (lu c <2K). The different T values observed by 
the two groups indicate that the speed of the oscilla- 
tory structure disappearing with temperature is sample 
dependent 4*^ To explain the temperature dependence 
the formation of an energy gap around the Fermi surface 
in the spectrum is suggested under microwave irradiation 
around the resistance minimal 

Our explanation of the temperature dependence of the 
magnetoresistance oscillations is based on the tempera- 
ture variation of the Landau level broadening T as deter- 
mined by Eq. 125|) . In a GaAs-based system, when the 
lattice temperature increases from around T = 1 K, the 
numbers of transverse and longitudinal acoustic phonons 
and thus the electron-phonon scattering strengths in- 
crease rapidly. In Fig. 7b we plot the zero-magnetic- 
field linear mobility due to transverse acoustic 
phonon scattering, /j,^ due to longitudinal phonon scat- 
tering, and the total mobility /iq as functions of lattice 
temperature T for the GaAs-based heterosystem with 
fio = 2.4 x 10 7 cm 2 V" 1 s" 1 at T = IK. We see that 
when temperature T rises from IK to 3 K the phonon 
related mobilities decline about an order of magnitude, 
leading to the total mobility fJ<o(T) decreasing about a 
factor of 2.2 and, according to l|25l) . T increasing about 
a factor of 1.5 (assuming a unchanged). The tem- 
perature growth of the Landau-level width due to this 
enhanced phonon scattering results in the strong tem- 
perature variation of the radiation-induced magnetore- 
sistance oscillation. The effect of such a mobility de- 
crease or the Landau level broadening on R xx is signifi- 
cant. Fig. 7a shows the calculated linear resistivity R xx 
as a function of lo/uj c at different lattice temperatures 
T = 1.0,1.5,2.0,2.5,3.0,3.5,4.0 and 4.5 K for the sys- 
tem of /zo(lK) = 2.4 x 10 7 cm 2 V~ 1 s _1 under a fixed mi- 
crowave illumination of frequency tojlix — 0.1 THz and 
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FIG. 7: (a) The longitudinal magnetoresistivity R xx at dif- 
ferent lattice temperatures T = 1.0,1.5,2.0,2.5,3.0,3.5,4.0 
and 4.5 K for a GaAs-based 2DEG subjected to a HF field 
E s sin(tirf) of frequency uj/2-k — 0.1 THz and amplitude 
E 3 — 80 V/cm. The system parameters are: electron den- 
sity N c = 3.0 x 10 n cm" 2 , dc mobility no = 2.4 x 10 7 cm 2 
V -1 s _1 at T = IK, and the broadening coefficient a — 12 
for all the curves, (b) Zero magnetic field linear mobility 
induced by transverse acoustic phonon scattering, /Uq , by 
longitudinal acoustic phonon scattering, fXg, and the total 
mobility fxg, for the system with /io(lK) = 2.4 x 10 7 cm 2 V _1 
s" 1 (dots), and for system with /Uo(lK) = 1.46 x 10 7 cm 2 V _1 
s _1 (solid line), (c) The amplitudes of the resistivity oscil- 
lation R xx (T) — R xx (5 K) at maxima and at minima versus 
temperature T for j = 1,2 and 3. (d) R XX (T)/ R XX (5K) is 
shown against 1/T on logarithmic scale for j = 2, 3 and 4. 



amplitude E s = 80 V/cm. The broadening parameter is 
fixed to be a — 12 for all the curves of different lattice 
temperatures, and the electron temperature is taken to 
equal to the lattice temperature in all the calculations 
in this section. The sensitive temperature dependence 
of the resistance oscillation is quite obvious. The mag- 
nitude of peaks and valleys of the oscillation diminish 
straightforward with increasing temperature from IK. 
At T > 4.5 K the oscillation structure almost disappear 
and R xx exhibits quite a flat form with changing lu/lo c 
for 7 C > 1.6. 

Comparing Fig. 7a with Fig. la^ we see that the role 
of decreasing inverse temperature at fixed radiation in- 
tensity is somewhat similar to that of reducing microwave 
intensity at fixed temperature. This similarity, how- 
ever, holds only within the range of T <4-5 K when di- 
rect phonon contribution to R xx is still less important 
than that of impurities. At higher temperature the di- 
rect phonon contribution may become dominant and the 
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FIG. 8: The longitudinal magnetoresistivity R xx at differ- 
ent lattice temperature T = 1.0,1.5,2.0,2.5,3.0,3.5,4.0,5.0 
and 6.0 K for a GaAs-based 2DEG subjected to a HF field 
E s sin(uit) of frequency lu/2tt = 0.1 THz and amplitude E s = 
70 V/cm. The system parameters are: electron density N e = 
3.0 x 10 n cm- 2 , dc mobility /i (lK) = 1.46 x 10 7 cm 2 V" 1 
s _1 , and the broadening coefficient a = 7.3 for all the curves. 
Also shown is the dark resistivity (E s = 0) at T — 6.5 K. 



structureless R xx will continue to grow with rising T and 
the R xx -~vs-j c curves of different T will not cross at inte- 
ger points 7 C = 2, 3, 4, 5, a feature of impurity scattering 
dominance. 

To show the temperature variation of R xx , we plot 
in Fig. 7c R xx — R XX (5K) at maxima and at minima of 
j = 1,2 and 3, as function of T. The oscillation dis- 
appears faster at larger j (lower magnetic field) than at 
smaller j (higher magnetic held). The positive values of 
Rxx/ Rxx(5 K) at peaks and at valleys also plotted against 
1/T on logarithmic scale. If we roughly ht the data with 
the form R XX (T) oc exp(-T /T), we have T » 10 K for 
j = 2, T w 4.5 K for j = 3, and T » 2.1 K for j = 4 in 
average over the range shown. 

The sensitivity of the temperature variation of the 
radiation-induced magnetoresistance oscillation is sample 
dependent. In GaAs-based systems the electron-phonon 
scattering strengths, thus the acoustic-phonon induced 
mobilities /i pt ' ) and /i P '' ) and their temperature behav- 
ior are essentially the same. Therefore, the tempera- 
ture variation of the total mobility /io depends mainly 
on the strength of the impurity scattering, which is al- 
most temperature independent within this range of T. 
The lowest curve in hg. 7b shows the T-dependence of 
the total mobility /xo (T) for the sample having T = 1 K 
mobility jUq(IK) = 1.46 x 10 7 cm 2 V _1 s _1 , which appar- 
ently exhibits a slower temperature change than that of 
Mo(lK) = 2.4 x 10 7 cm 2 V -1 s -1 system. 

In Fig. 8 we illustrate the linear resistivity R xx as a 
function of lu/lu c at different lattice temperatures T = 
1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 5.0 and 6.0 K for the system 
of Mo(lK) = 1.46 x 10 7 cm 2 V -1 s _1 under a hxed mi- 
crowave irradiation of frequency w/2tt = 0.1 THz and 
amplitude E s — 70 V/cm. The broadening parameter 
is fixed to be a — 7.3 for all the lattice temperatures. 



The speed of the oscillatory structure disappearing with 
rising temperature is apparently slower than the system 
shown in Fig. 7a. We see that at T > 6.0 K the oscillation 
structure essentially disappear and R xx approaches a flat 
form of dark (E s = 0) curve of corresponding tempera- 
ture (see Fig. 8) for j c > 1.6. Note that different from 
T = 1 K case, where impurity scattering dominates, at 
T = 6K the acoustic phonon scatterings already yield 
a substantial contribution to both zero-field mobility fio 
and to radiation-induced magnetoresistivity R xx - 



VI. CONCLUSION 

Based on the balance-equation model for magneto- 
transport in Faraday geometry, we have carried out a 
detailed theoretical investigation on microwave-radiation 
induced magnetoresistance oscillations recently discov- 
ered in high-mobility GaAs-based two-dimensional elec- 
tron systems. We find that for systems having zero- 
field linear mobility /j (1K) < 2.4 x 10 7 cm 2 V^ 1 s _1 
multiphoton-assisted impurity scatterings are the main 
mechanisms responsible for radiation-induced magne- 
toresistance oscillations at temperature T <4K. The am- 
plitude of the R xx oscillation grows roughly following 
the microwave power under weak illumination and fol- 
lowing the microwave amplitude under medium illumi- 
nation before it saturates even decreases with continu- 
ing increase of the microwave strength under strong ir- 
radiation. It is shown that the strongest oscillations ap- 
pear in the linear longitudinal magnetoresistance and a 
finite dc current bias always suppresses the oscillation. 
We find that, different from the SdH oscillation which 
is easily suppressed by a few-degree rise of the elec- 
tron temperature, the radiation-induced magnetoresis- 
tance oscillations are quite insensitive to the modest elec- 
tron heating as long as the lattice temperature remain the 
same. Although the magnetoresistivities directly stem- 
ming from photon-assisted transverse and longitudinal 
acoustic phonon scatterings also exhibit pronounced os- 
cillations under microwave irradiation, they contribute 
only a small part of the total R xx in the temperature 
range of T < 4 K. Nevertheless, as we proposed, that it is 
just this acoustic phonon scattering that gives rise to the 
sensitive lattice temperature dependence of radiation- 
induced resistance oscillations right from T = 1 K. We 
showed that the growth of the Landau level broadening 
resulting from the enhancement of acoustic phonon scat- 
terings with increasing lattice temperature leads to the 
observed temperature suppression of the oscillation. 
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